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ABSTRACT 


The  symmetry- plane  laminar  boundary  layer  of-  an  impulsively-started 
ellipsoid  of  revolution  at  high  incidence  is  solved  to  shed  light  on  some 
basic  characteristics  of  three  dimensional,  unsteady  flows.  The  governing 
equations  are  formally  similar  to  those  for  the  three-dimensional,  steady  case, 
so  the  same  method  of  solution  and  computer  programs  previously  developed  were 
employed  in  the  present  work.  The  most  important  result  obtained  is  concerned  • 
with  the  meridional  skin  friction  cf^  .  The  zero-c^  point  (at  which  cf^  vanishes  ) 
does  not,  as  expected,  move  forward  as  time  increases,  instead  it  remains  over  the 
rear  body.  As  t-*» 00  ,  it  jumps  to  the  front  nose.  This  implies  that  there  is  no 
flow  separation  over  the  symmetry-plane  at  finite  times.  In  the  meanwhile, . it 
is  argued  that  separation  must  occur  on  two  sides  of  the  body.  This  situation 
leads  us  to  propose  a  new  unsteady  separation  sequence,  i.e.  an  open  type  separation 
prevails  at  earlier  times,  while  a  closed  type  of  separation  occurs  only  at  the 
steady-state  condition.  This  sequence  presents  a  sharp  contrast  to  the  conventional 
notion  of  unsteady  sfipaaraEtion  which  consists  of  a  series  of  closed  separations  only. 
Furthermore,  this  sequence  with  the  time  as  the  parameter  is  found  to  be  similar 
to  that  for  previously-studied  steady  flows  with  varying  incidences. 
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1 .  INTRODUCTION 


The  complexity  of  boundary  layer  research  increases  from  two-dimensional 
steady  problems  to  two-dimensional  unsteady,  three-dimensional  steady  and  finally 
to  three-dimensional  unsteady  problems .  Considerable  progress  has  been  made 
in  recent  years  to  two-dimensional  unsteady  and  three  dimensional  steady  problems, 
but  the  three-dimensional,  unsteady  area  remains  to  be  the  last  frontier  to  be 
fully  exploited.  Three-dimensional,  unsteady  solutions  are  very  scarce.  A 
classical  work  in  this  subject  is  due  to  Squire  (1954)  who,  using  a  time-series 
expansion  method  investigated  the  motion  of  an  impulsively-started  ellipsoid. 

Recent  publications  in  this  area  include,  for  example,  those  of  Dwyer  and  Sherman 
(1979) and  Williams  (1980). 

As  a  prelude  to  full  three-dimensional  investigations,  the  present  work 
considers  the  symmetry-plane  boundary  layer  over  an  ellipsoid  of  revolution.  The 
steady  counterpart  of  this  problem  was  studied  before  by  one  of  the  present  authors 
(Wang  1970)  and  later  by  Cebei  et  al  (1980) .  Although  consideration  was  restricted 
to  the  symmetry-plane,  the  results  revealed  features  of  fundamental  significance 
for  (steady)  three-dimensional  flows  in  general  (Wang  1974a,  1976).  It  was  our 
original  motivation  that  the  present  work  will  achieve  the  same  for  the  corresponding 
unsteady  case.  Our  results  to  be  presented  later  seem  to  bear  out  this  expectation. 

The  specific  unsteady  motion  considered  is  that  of  an  impulsive  start.  This 
problem  provides  an  interesting  contrast  to  the  classical  impulsively-started 
circular  cylinder  problem.  In  the  latter  case,  the  zero-skin-friction  point  starts 
from  the  rear  stagnation  point,  moves  gradually  forward  as  time  increases  and 
approaches  finally  the  steady-state  separation  position  (0  *  104.5°).  For  the 
♦  present  ellipsoid  case,  the  zero-skin-friction  point  on  the  leeslde  symmetry-plane 

was  at  first  thought  to  move  very  much  the  same  way.  From  our  previous  investiga¬ 
tions  (Wang  1974b) ,  the  steady-state  position  bf  the  zero-skin-friction  point  on 
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the  symmetry-plane  is  located  at  the  front  nose,  while  a  closed  type  of  separa¬ 
tion  prevails  over  the  whole  body  (figure  lc) .  Consequently  if  the  zero-skin- 
friction  point  on  the  leeside  symmetry-plane  gradually  moves  forward  as  expected  , 
then  it  seems  natural  to  think  that  similar  overall  separation  pattern  would  also 
prevail  over  the  body  at  earlier  times  as  well .  This  leads  to  an  unsteady  sequence 
of  developments  as  depicted  in  figures  l(a)-(c),  each  figure  displays  a  closed  type 
of  separation  except  the  separated  area'  expands  with  time.  This  sequence  represents, 

i 

in  fact,  the  conventional  idea  of  unsteady  separation.  However,  in  Section  4, ‘we 
shall  show  that  our  present  results  did  not  turn  out  the  way  as  expected  and  in 
Section  5,  we  shall  propose  instead  a  new  sequence  of  unsteady  separation  consisting 
of  both  open  and  closed  types  of  separation  patterns. 

Finally  we  would  like  to  comment  on  the  current  status  of  the  boundary  layer 
theory  relative  to  the  "thin  layer"  theory  (Pulliam  and  Stegar  1980,  Hsieh  1981). 

The  basic  idea  of  the  "thin  layer"  approximation  is  similar  to  that  of  the  boundary 
layer  theory  except  the  pressure  is  not  imposed  according  to  an  inviscid  solution, 
instead  it  is  calculated  along  with  the  boundary  layer.  This  allows  the  inviscid 
boundary  layer  interaction  and  hence  broadens  the  range  of  validity  beyond  that  of 
the  boundary  layer  theory.  However  a  thin-layer  problem  is  mathmatically  an 
elliptical  problem  in  the  space  coordinates,  whereas  a  boundary  layer  problem  is 
parabolic.  Solving  an  elliptical  problem  requires  larger  computer  storage  and 
longer  computation  time'  so  that  adequate  resolutions  were  found  not  possible  on 
a  CDC  7600  computer. 

This  difficulty  led  later  to  the  idea  of  parabolicization,  i.e.  treating 
an  elliptic  problem  in  a  parabolic  fashion  (Schiff  and  Steger  1980).  This 
procedure  requires  that  there  is  no  upstream  pressure  influence  in  the  main  flow 
direction.  The  latter  condition  is  best  satisfied  in  a  supersonic  cone  flow  for 
which  the  longitudinal  (or  meridional)  perssure  gradient  is  identically  zero. 
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However  even  in  this  particular  case,  the  above  condition  is  not  completely  met 
because  the  flow  is  no  longer  strictly  conical  when  the  inviscid-boundary  layer 
interaction  is  considered.  In  a  general  supersonic  flow,  there  is  always  a. 
subsonic  part  of  the  boundary  layer  where  the  above  condition  cannot  be  satisfied. 

For  subsonic  flows,  parabolicization  cannot  be  justified.  Thus  thin-layer  theory 
is  more  general,  but  is  accompanied  by  the  above  shortcomings.  At  the  present 
state  of  the  art,  it  appears  that  boundary  layer  theory  is  still  the  most  appropriate 
approach  for  general  three-dimensional  viscous  flows,  because  even  to  this  lower 
approximation  (compared  to  the  thin-layer  version) ,  unsteady  solutions  are  still  not 
available.  On  the  other  hand,  one  may  favor  the  thin-layer  version  when  computer's 
storage  and  time  are  not  of  much  concern  under  certain  circumstances,  because  it 
avoids  to  solve  the  inviscid  and  viscous  flows  separately  aside  from  the  interaction 
between  them  mentioned  above. 

I 
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2.  EQUATIONS 

The  governing  equation  for'  the  present  symmetry-plane  boundary  layer 
written  in  non-dimensional  form  consists  of 


Bu 

4. 

Bu 

4. 

Bu 

.  -BP 

-+  iii 

Bt 

r 

Uh  By 

P 

T* 

w  b7  ’ 

h  Bu 

u 

3z2 

* 

3v9 

4- 

3V0 

x. 

w  3Vq 

2 

+ 

•+  ^ 

dr 

St 

T 

uhyBu 

T 

r 

r 

hy3u 

Bu 

X  1 

u  3r 

4. 

Bw 

—  n 

575 

T 

rhy3u 

r 

Bz 

*  u . 

rS  9 


3  v9 

.  2, 
3  z 


(la) 

(lb) 


(lc) 


Aside  from  the  unsteady  terms ,  these  equations  are  otherwise  identical  to  these 
for  the  steady  case  treated  before  (Wang,  1970).  The  same  notations  are  also  used. 
Briefly,  referring  to  .figure  2, a  is  the  incidence  angle,  V  and  9  are  two  surface 
coordinates,  z  is  the  normal  coordinate,  u,  v  and  w  are  the  corresponding  velocities, 
U  and  V  are  the  inviscid  velocities  at  the  outer  edge.  hu  and  r  are  the  metric 
coefficients  and  r,  in  fact,  is  also  the  cross-sectional  radius. 

hy  =  [(1  -  e2  u2)  /  (1  -  u2)f, 

r  =*  [(1  -  e2)  (1  -  u2)  ]% 

where 

2  2  ^ 

e  *  (1  -  b  /a  ) 

a  and  b  are  the  semi-major  and -minor  axis  of  the  ellipsoid. 

At  the  plane  of  symmetry,  v  3  0,  so  that  v0(=  9v/B9)  is  taken  as  an 
independent  variable  along  with  u  and  w.  The  boundary  conditions  are 


«  ■  D,  v  *  V  (i.e.  BV/B0)  at  z-*-°°  , 

u  -  vQ  -  0  at  z  -  0, 

where  U,  V  and  the  pressure  gradients  are  known; 
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(2a,b,c) 
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Based  on  the 
defined  by 


u-velocity,  the  skin  friction  and  the  displacement  thickness  are 
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Analogously,  we  define  based  on  the  -  profile, 


(c  )  - 

V  f  V0  fg-  \  dz  )  z+0  , 

Ave*  Jf(1-Vve)dz  • 


(4c, d) 


Although  the  vQ-  profile  has  always  been  calculated  in  the  symmetry-plane 

* 

type  of  boundary  layer  studies,  analogous  displacement  thickness  A  ^  was, 
however,  not  considered  except  by  Wang  (1971a).  Even  there  A  ^  was  not  recog¬ 
nized  to  be  needed  in  evaluating  the  total  displacement  thickness. 

Moore  and  Ostrach  (1957)  derived  the  equation  for  the  unsteady  total 

A* 

displacement  thickness  “  , 
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In  applying  this  equation  to  the  symmetry- plane  case,  care  must  be  taken  to 


incorporate  Avg  into  the  final  expression; 
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In  the  literature,  only  has  been  reported.  Here  we  shall  present  results 
of  both  A^  and  Av^,  but  the  calculation  of  A  is  still  not  carried  out  yet. 


3.  METHOD  OF  SOLUTION 


The  equations  (la,b,c)  are  similar  in  structure  to  those  for  three- 
dimensional,  steady  boundary  layers  investigated  before  (Wang  1974a).  So  the 
same  numerical  methods  and  computer  programs  developed  earlier  were  used  for 
the  present  work. 

3.1.  Initial  Profiles. 

To  start  the  computation,  two  sets  of  initial  profiles,  temporal  and  spatial, 
are  required. 

Temporal  Initial  Profiles  (t=tQ).  Immediately  after  (impulsive)  starting, 
potential  flow  prevails.  For  the  temporal  initial  profiles , the  usual  first  approxi¬ 
mation  of  a  small-time  expansion  is  adequate,  i.e. 


(6a) 

(6b) 


where  erf  is  the  error  function. 


Spatial  Initial  Profiles  (y=yQ).  Based  on  the  argument  that  diffusion  and 
convection  affect  the  downstream  flow  only,  so  that  solutions  near  the  stagnation 
point  are  nearly  independent  of  time.  Thus  the  required  spatial  initial  profiles 
may  be  provided  by  the  steady  symmetry-plane  solutions  near  the  stagnation  point. 


3.2.  Dependence  Zone 

For  steady  three-dimensional  boundary  layers,  calculation  was  found  to  have 
to  follow  the  zone  of  dependence  rule  (Raetz  1957,  Wang  1971b).  For  steady,  symmetry- 


plane  problems,  the  wedge-shaped  dependence  zone  shrinks  to  the  symmetry-plane  itself 
and  it  was  found  (Wang  1970)  that  calculation  stops  at  the  onset  of  the  reversal  of 
the  u-velocity  profile,  but  reversal  of  the  v^-profile  was  unexpectedly  calculated 
with  no  sign  of  any  difficulty.  This  result  was  at  first  looked  upon  with  suspicion, 
but  later  calculation  of  the  full  three-dimension  problems  as  well  as  the  repeat 
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calculation, of  the  same  symmetry-plane  problem  confirmed  that  result  to  be 
correct.  However,  a  theoretical  explanation  was  not  available. 

For  unsteady  cases,  two-dimensional  or  the  symmetry-plane ,  an  analogous 
zone  of  dependence  (Wang  1975,  1979)  must  be  observed.  So  at  the  beginning  of 
the  present  work,  we  were  faced  again  with  the  same  question;  what  are  precisely 
the  differences  in  logic  for  the  calculation  of  the  u-profile  and  the  v  -profile? 
In  the  unsteady  case,  we,  in  fact,  further  know  that  reversal  of  u-profile  can  be 
calculated  up  to  separation  unlike  in  the  steady  case  where  calculation  of  the  . 
reversal  of  the  u-profile  is  prohibited  (Wang  1979) . 

To  provide  the  missing  explanation,  we  carried  out  a  characteristics  study 
of  the  sys.temof  equations  (la,b,c).  The  determinants  (Wang  1971b)  of  the  character¬ 
istics  and  subcharacteristics  were  found  respectively  to  be 


n  _  2  /  \ 


_  3S1/3!!,  u  d9.\2 

Q  \TE  +  v  9U  +  w  j  • 


(7a) 


(7b) 


which  shows  that  the  addition  of  equation  (lb)  for  v  (compared  to  equations  for 
two-dimensional  cases)  merely  increase  the  power  from  3  (Wang  1975)  to  5  in  equation 
(7a)  and  from  1  to  2  in  equation  (7b).  Otherwise  the  same  dependence  rule  applies 
to  both  the  two-dimensional  and  the  symmetry-plane  problems.  The  underlying  reason 
why  reversal  of  v  is  not  prohibited  by  the  dependence  rule  is  recognized  due  to 

v? 

the  fact  that  in  both  equations  (la,b),  the  differential  operator 


3  4.  X  9 

3t  u  lTliZ  +  w  3z 

does  not  involve  v^.  Yet  it  is  this  operator  which  determines  the  subcharacteristics 
in  equation  (7b). 

3.3.  Computational  Details 

The  initial  time  t  was  set  to  be  0.01  and  the  initial  space  U  was  taken 

n  D 


J  4. 


to  be  -0.984,  very  close  to  the  front  vertex  (u  =  -1.0).  With  initial  profiles, 
calculation  marches  from  the  front  stagnation  point  toward  the  rear  body  at  a 
fixed  time  t.  At  each  station  (u,t),  solution  between  the  body  (z=0)  and  the 
outer  edge  (z-*-  °°)  is  implicitly  calculated.  Then  the  same  process  is  repeated 
as  time  increases  from  t  to  t+  At  until  the  whole  problem  is  completed. 

Ar  high  incidaice  (45°),  large  pressure  gradients  prevail  in  both  the  meri¬ 
dional  and  circumferential  directions,  especially  over  the  two  ends.  Computational 
difficulties  due  to  poor  convergence  were  experienced  at  the  beginning  of  the  ' 
present  work,  it  was  only  diagnosed  later  that  unusually  small  steps  must  be  employed 
Figure  3  compares  the  results  of  two  calculations  with  different  mesh  sizes.  Although 
At  =  Au®  0 . 005  is  normally  considered  to  be  small,  results  still  show  large  deviations 
Immediately  after  the  impulsive  start  (t<0.05),  relatively  large  meshes 
(At  =  0.005,  Au  =  0.005)  could  be  used.  Then  smaller  space  step  (Au  ■  0.00125  was 
required  near  two  ends  of  the  body,  finally  At  must  also  be  reduced  (0.0025  and 
0. 00125). For  0.05£ t<_ 0.16,  our  computation  covers  the  whole  leeside  symmetry-plane; 

At  =  0.00125,  but  two  different  Ay ' s  were  used  to  save  the  computing  time.  These 
values  are  0.00125  fcr  two  ends,  and  0.005  for  the  mid-body.  For  0.16<t<1.0,  only 
the  front  nose  portion  was  calculated  (reason  is  given  in  section  4),  At  ■  0.00125, 

Au  =  0.00125  and  even  0.000625.  Such  three-dimensional  (t,u  z)  computations  with 
small  steps  are  very  expensive.  Compared  to  that  for  the  u-profile  ,  calculation 

of  the  v  -profile  especially  demands  finer  steps  in  order  to  satisfy  the  outer- 
y 

edge  condition  (i.e.  to  merge  smoothly  with  V  ).  This  is  attributed  to  the  lateral 

v? 

(circumferential)  pressure  gradient  being  more  severe  than  the  meridional  pressure 
gradient  for  the  problem  under  consideration. 
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4.  RESULTS  OF  HIGH  INCIDENCE 

The  windside  symmetry-plane  problem  is  more  straightforward,  so  only  the 
leeside  results  (cx«  45°)  are  presented  in  details.  For  low  to  moderate  incidences, 
the  unsteady  growth  is  expected  to  have  a  very  different  structure,  and  will  be 
reported  separately  later  on. 

4.1  Pressure  Variation 

Figure  4  gives  the  variation  of  pressure  p,  longitudinal  pressure  gradient 

2  2 

9p/h^9y  and  lateral  pressure'  curvature  s p/r9  9.  The  last  two  are  directly 
involved  in  equations  (lb,c).  The  pressure  gradient  is  fairly  mild  over  the  mid¬ 
body,  but  extremely  large  near  two  ends.  This  suggests  that  computational  difficulties 
as  well  as  the  most  interesting  features  of  fluid  mechanics  would  likely  occur  near 
the  ends  also. 

4.2  Skin  Friction  c£ 

f  U 

Figure  5(a)  gives  the  meridional  skin-friction,  cf ^  ,  along  the  leeside 
symmetry-plane.  At  first,  cf  decreases  from  the  front  to  become  zero  at  the 
rear  stagnation  point.  But  as  time  increases,  two  important  features  emerge: 

(1)  The  zero-Cj  po 'nt  does  not  gradually  move  forward,  instead  it  remains  over 
the  rear  body;  (2)  there  develops  a  minimum  point  at  the  front,  and  furthermore, 
this  minimum-cF  continues  to  decrease  until  finally  it  becomes  zero.  At  that 
time,  the  zero-c^  point  suddenly  jumps  from  the  rear  to  the  front. 

Calculationsover  the  entire  leeside  symmetry-plane  were  carried  out  up  to 
t  *  0.16.  At  that  stage,  it  was  felt  that  the  abovementioned  two  features  were 
established,  and  the  remaifting  task  was  only  to  determine  how  fast  this  minimum- 
Cj^  will  approach  zero.  For  this  purpose,  it  is  unnecessary  to  calculate  over 
the  entire  symmetry-plane.  Hence  later  calculations  were  limited-  to  the  front 
body  only  in  order  to  save  the  computing  cost.  The  details  of  c^  from  the  later 
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calculations  are  shown  in  figure  5(b).  As  time  increases,  it  is  seen  that 

approaches  zero  at  an  extremely  slow  rate.  For  example,  between  t  =  0.5  and 

1.0,  (c,  )  is  onlv  reduced  from  0,2403  to  0.1828.  With  At  =  0.0025,  this 

f  U  mm 

part  of  the  calculation  alone  took  200  cycles  with  £p  =  0.000625.  The  slow  change 
led  us  to  interpret  that  (cfy)m^n'  approaches  zero  asymptotically  as  t-*«»  . 

4.3.  Analogous  Skin  Friction  (cf)VQ 

Figures  6 (a) -(b)  give  the  corresponding  results  of  the  analogous  skin  friction 
(c^)^  based  on  the  profile  of  vQ.  (c^)^- determines  among  others  the  sign  of. 
v-velocity  near  but  off  the  symmetry-plane.  Change  of  the  sign  of  (c^)  ^  means 
the  reversal  of  the  v-velocity .  At  small  times  (t<0.10),  figure  6(a)  shows  the 
point  of  zero-(cpvg,  Pi,  moving  gradually  forward  from  the  rear  stagnation 
point.  At  0.10<t<0.13,  additional  v-reversals  occur  at  and  over  the  front. 

As  time  further  increases,  R2  moves  rearward  and  merges  with  R^,  then  disappearing 
altogether.  At  t>0.13,  there  exists  only  R^  at  the  front.  The  above  description 
of  the  zero-(c^)vQ  point  is  summarized  in  figure  6(c)  which  displays  clearly  its 
non-mono tonic  behavior.  At  0.10£t<0.13,  there  are  R^,  R2  and  R^ (corresponding  to 
t=  0.12).  As  time  increases ,  R^  moves  from  D  (i.e.  rear  stagnation  point)  to  C,  R9 
from  B  to  C  and  R^  from  B  toward  A.  The  curve  ABCD  divides  the  u,t-plane  into 
two  parts;  i.e.  whether  v-velocity  near  the  symmetry-plane  is  reversed  or  not. 

The  surface  flow  pattern  near  the  symmetry-plane  in  accordance  with  figure  5 
and  6  are  illustrated  in  figure  7.  In  figure  7(a),  the  flow  at  earlier  times 
points  back  and  up  (except  between  R^  and 

tion  point),  whereas  in  figure  7(d),  the  flow  points  back  and  down.  Figures 

7(b) -(c)  indicate  the  intermediate  changes.  Figure  7(e)  suggests  that  after 

the  point  of  zero-c^  jumps  to  the  front,  the  steady-state  separation  point 

S  (U«  -0.921)  is  located  downstream  of  R,  (u*  -0.955).  This  makes  S  as  a 
s  j  s 

saddle  point  of  separation  as  expected.  More  about  separation  will  be  discussed 


A^  where  A^  stands  for  the  rear  stagna- 


in  section  5.  In  any  case,  the  changes  of  v-velocity's  sign  back  and  forth  is 
rather  unusual,  indicating  the  complex  structure  involved  in  an  unsteady  development. 

4.4  B, sic  Profiles 

The  reasons  about  why  c,  and  (c,,)  _  show  unusual  variations  can  all  be 

t  u  f  v0 

traced  out  from  the  results  of  the  corresponding  basic  profiles.  Figure  8(a) 
shows  the  u-velocity  profiles  at  t  =*  0.12.  The  non-monotonic  variation  of  cf 

E  y 

shown  in  figure  5(a)  is  a  reflection  of  the  intersection  of  the  profiles  at 
different  u~stations.  Near  the  body  (z<0.4),  the  curve  f  or  y  *  0.716  stands 

between  those  forU*  -0.946  and  -0.133.  This  implies  that  c,  is  the  smallest 

fy 

at  the  front  body  y  =  -0.946,  increases  at  the  mid-body  U=  -0.133  and  decreases 
at  the  rear  body  0.716.  Figure  8(b)  shows  the  profile  of  vQ  at  t  *  0.12. 

That  the  curve  for  u»  -0.133  falls  between  those  for  y»  -0.946  and  -0.833  (rather 
than  between  U»  -0.833  and  0.716)  is  again  responsible  for  the  sign  of  v-velocity 
to  change  back  and  forth . 

4.5  Displacement  Thickness 

Figure  9  shows  the  displacement  thicknesses  A^  and  A*Q  .  Unlike  in  the 
normal  situation,  these  thicknesses  do  not  grow  montonically  in  the  downstream 
direction,  instead  they  increase  sharply  near  the  front  nose  and  then  curve 


downward  immediately  thereafter.  This  phenomenon  again  follows  from  the 


profiles  shown  in  Figures  8(a)-(b),  and  was  similarlv  noted  before  in  the 

steady  cases  (Wang  1970,  1971a).  A*  increases  much  faster  than  A*  with  increas- 

vg  U 

ing  time,  reflecting  the  dominance  if  the  circumferential  flow.  A*0  is,  in  fact, 
presented  here  for  the  first  time,  its  importance  has  been  to  date  overlooked  in 


the  literature. 
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5.  IMPLICATION  OF  3-D  UNSTEADY  SEPARATION 
The  variation  of  the  meridional  skin-friction  c^  discussed  in  connection 
with  figure  5(a)  has  an  important  bearing  with  respect  to  three-dimensional 
separation , i . e .  separation  on  the  whole  body  rather  than  along  the  symmetry-plane 
itself.  A  similar  situation  was  noted  previously  in  the  steady  case  (Wang  1971a) 
from  which  figure  10  is  reproduced  to  show  the  corresponding  skin  friction  c^. 

As  the  incidence  increases,  there  develops  a  minimum  point  of  c^  over  the  front 
nose.  This  minimum  dips  deeper  until  it  becomes  zero.  When  that  happens,  the 
zero-skin-friction  point  jumps  from  the  rear  bodv  to  the  front  nose. 
Comparing  figures  10  and  5(a),  it  is  clearly  seen  that  in  both  cases,  the  zero- 
skin-friction  point  remains  over  the  rear  body,  while  a  minimum  develops  over  the 
front  nose.  The  varying  parameters  are  the  time  in  figure  5(a)  and  the  incidence 
in  figure  ]0  . 


The  behavior  of  c^  in  figure  10  led  us  later  to  conceive  the  open  vs  closed 
separation  concept  (Wang  1972,  1974a,  1976).  The  basic  idea  of  a  closed  separation 


fits  well  with  the  traditional  notion  of  separation,  but  that  of  an  open  separation 
proves  to  be  just  the  opposite.  In  the  past  few  years,  the  open  separation  idea 
has  been  confirmed  by  many  experiments  and  numerical  solutions  (Wang  1976,  1982). 

By  the  same  reasoning,  we  propose  here  an  unsteady  counterpart  of  the  open-vs- 
closed  separation  sequence.  The  arguments  are  as  follows: 

For  the  steady  two-dimensional  case  or  the  steady  symmetry-plane  case 
(Wang  1970) ,  separation  is  marked  by  the  vanishing  of  the  meridional  skin  friction 
c^  .  For  the  unsteady,  two-dimensional  case,  it  was  demonstrated  (Wang  1979)  that 
separation  is  identified  by  the  running- together  of  the  analogous  limiting  stream¬ 
lines  in  the  U,  t-plane,  the  same  idea  holds  for  the  unsteady,  symmetry-plane 
problem  concerned  here.  By  this  criterion,  although  the  vanishing  of  c^y  is  not 
synonymous  to  separation,  yet  separation  occurs  only  after  c^  goes  to  zero  and 


14 

then  becomes  negative.  However,  in  section  4,  since  we  noted  already  that 
the  zero-c.^  point  does  not  continue  to  move  forward,  so  we  conclude  that 
there  is  no  flow  separation  on  the  leeside  symmetry-plane.  Furthermore,  in 
order  to  have  the  separation  patterns  of  figure  lb,c,  there  must  be  a  zero- 
point  on  the  symmetry-plane.  This  condition  is  again  denied  by  the 
present  results,  so  we  conclude  that  the  sequence  shown  in  figure  l(a)-(c) 
is  not  to  happen. 

We  must  next  emphasize  that  lack  of  separation  on  the  symmetry-plane 
does  not  mean  there  is  no  separation  on  the  sides  of  the  body.  In  fact,  due 
to  the  fact  that  the  circumferential  flow  over  most  of  the  body  is  subject  to 
larger  pressure  gradient  than  the  meridional  flow,  so  that  a  mainly  lateral 
separation  must  occur  on  two  sides  of  the  symmetry-plane  at  much  earlier  times. 
Besides,  we  knew  beforehand  that  the  steady-state  separation  at  such  a  high 
incidence  is  of  a  closed  type  as  shown  in  figures  lc  or  11c.  It  is  inconceivable 
that  such  a  whole  closed  separation  pattern  is  established  instantaneously  at  the 
moment  when  the  zero-c^ u  point  jumps  to  the  fore  body.  Based  on  these  reasons, 
we  propose  that  the  likely  sequence  of  developments  is  as  depicted  in  figures 
lla,fc,c.  In  other  words,  an  open  type  of  separation  starts  first  (t=t^)  over 
the  rear  body,  extends  forward  as  time  increases  (t=t2)  until  the  open  separation 
line  finally  intercepts  the  symmetry-plane  upon  the  jump  of  the  zero-c^  point 
from  the  rear  body  to  the  front,  thus  completing  the  formation  of  a  closed  separation 
for  the  steadv-state  condition. 

The- proposed  unsteady  sequence  of  separation  is  analogous  to  a  steady  counter¬ 
part  suggested  by  one  of  the  present  authors  (Wang  1976) .  While  the  latter  has 
been  confirmed  by  experiments  and  calculations,  we. expect  that  experimental 
demonstration  of  the  unsteady  separation  will  be  more  difficult,  whereas  full 
three-dimensional  unsteady  boundary  layer  solutions  also  remain  to  be  reported. 
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The  preceding  conclusion  of  no  separation  of  finite  times  touches  on  an 
interesting  story.  Over  the  past  few  years,  there  was  a  controversy  (Wang  1981) 
on  this  question  in  the  literature.  Cebeci.  (1979)  advanced  a  similar  idea  in 
connection  with  an  impulsively-started  circular  cylinder  problem.  His  idea 
was  disputed  by  other  researchers  (including  one  of  the  present  authors,  Wang)  and 
has  since  been  known  to  be  in  error  for  the  cylinder  case.  But  ironically,  our 
present  results  do  seem  to  provide  an  example  in  agreement  with  his  idea  as  far 
as  the  separation  on  the  symmetry -plane  is  concerned..  The  circumstances  here. are 
of  course  different  from  what  he  originally  conceived,  the  symmetry-plane  separation 
is  only  a  part  of  the  overall  three-dimensional  separation,  whereas  his  problem  was 
a  two-dimensional  one. 
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6.  CONCLUSIONS 

For  an  impulsively-started  symmetry-plane  problem  at  high  incidence,  the 
circumferential  flow  near  the  symmetry  plane  was  found  to  reverse  its  direction 
back  and  forth. 

The  variation  of  the  meridional  skin  friction  c^  is  surprisingly  found  to 
be  similar  to  that  for  the  corresponding  steady  case  studied  before,  although 
the  varying  parameter  is  the  time  in  the  present  unsteady  case,  while  the  same 
is  the  incidence  angle  in  the  previous  steady  case.  The  zero-c^point  was 
found  to  remain  over  the  rear  body  at  all  finite  times  until  it  jumps  to  the 
front  nose  to  reach  the  steady-state  condition  at  t-*-». 

Since  the  zero-c^  point  does  not  gradually  move  forward,  it  is  concluded 
that  there  is  no  Reparation  on  the  symmetry-plane  at  all  finite  times.  Instead 
an  open  separation  is  likely  to  develop  first  over  the  sides  of  the  rear  body, 
and  gradually  extends  forward.  When  the  zero-cf^point  jumps  to  the  front  nose 
and  the  open  separation  line  intersects  the  symmetry-plane,  then  a  closed  separ¬ 
ation  is  formed  at  the  steady-state  condition.  On  the  one  hand,  this  open  and  closed 
separation  sequence  provides  a  striking  contrast  to  the  conventional  notion  of 
unsteady  sequence  consisting  of  closed  separations  only.  On  the  other  hand,  this 
sequence  with  time  as  the  parameter  is  found  to  be  similar  to  that  for  steady 
flows  with  varying  incidences.  Our  open  vs.  closed  separation  idea  originally 
conceived  for  steady  flows  is  now  extended  to  unsteady  flows. 
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Figure  3.  Comparison  of  different  meshes 
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